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Abstract The incompressible boundary layer in the axial flow past a cylinder has been shown Tutty et. 
al-dTi I) to be stabler than a two-dimensional boundary layer, with the helical mode being the least stable. 
In this paper the secondary instability of this flow is studied. The laminar flow is shown here to be always 
stable at high transverse curvatures to secondary disturbances, which, together with a similar observation 
for the linear modes implies that the flow past a thin cylinder is likely to remain laminar. The azimuthal 
wavenumber of the pair of least stable secondary modes (m+ and m_) are related to that of the primary (n) by 
m+ = 2n and m_ = ~n. The base flow is shown to be inviscidly stable at any curvature. 
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1 Introduction 

At low to moderate freestream disturbance levels, the first step in the process of transition to turbulence in a 
boundary layer is that at some streamwise location, the laminar flow becomes unstable to linear disturbances. 
While this instability and the events that follow have been investigated in great detail for two-dimensional 
boundary layers during the past century, much less work has been done on its axisymmetric counterpart, the 
incompressible boundary layer in the flow past a cylinder, notable exceptions being the early and approximate 
linear stability analysis of Rao|7| and the more recent and accurate one of Tutty et. al.| 1 1]. In Rao's work, the 
equations were not solved directly and the stability estimates had severe limitations. Tutty et. al.| 1 1 1 showed 
that non-axisymmetric modes are less stable than axisymmetric ones. The critical Reynolds number was found 
to be 1060 for the n = 1 mode and 12439 for n = 0. The instability is thus of a different character from that 
in two-dimensional boundary layers, since Squire's (1933) theorem, stating that the first instabilities are two- 
dimensional, is not applicable in this case. The expected next stage of the process of transition to turbulence, 
namely the secondary modes of instability, have not been studied before, to our knowledge, although the 
turbulent flow over a long thin cylinder has been studied by Tutty 1 10], who computed the meanflow quantities 
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Fig. 1 Schematic diagram siiowing tiie coordinate system. 



to within experimental accuracy. Practical applications, on the other hand, are numerous. For example, the 
axial extent of turbulent flow would determine the signature that submarines leave behind themselves, apart 
from the drag. The transition to turbulence over the bodies of large fish too would be partially controlled by 
transverse curvature. 

The secondary instability of incompressible laminar flow past an axisymmetric body is thus the focus of 
this paper. We present results only for the flow past a cylinder, but the equations derived here and the solution 
method can be used for arbitrary axisymmetric bodies. We show that the overall effect of transverse curvature 
on incompressible boundary layers is to stabilise secondary disturbances. Remarkably no instability is found 
at any Reynolds number at higher curvatures, i.e., when the boundary layer thickness becomes comparable to 
the body radius. This implies that the boundary layer past a thin cylinder would tend to remain laminar, or to 
relaminarise downstream even were it to go turbulent. 

We note in contrast that longitudinal curvature, and the resulting Gortler (1940) vortices on concave walls, 
have been well studied (see e.g. Benmalek and Saric UJ) and so have instabilities in three-dimensional bound- 
ary layers due to streamline curvature E). 



2 Mean flow 



The unperturbed laminar flow is obtained by solving the incompressible steady boundary layer equation for 
the axial component of velocity: 

^dU ^dU l/d^U ldU\ 
ox or R\ or^ r or J 

together with the continuity equation 
dU dV V ^ 

^ + ^ + -=0, (2) 
ax or r 

and the boundary conditions 

U{Q,r) = l, U{x,l)=Q, y(jc,l)=0 and [/(x,oo) = l. (3) 



Here the x coordinate is along the surface of the body and r is normal to the body surface and measured from 
the body axis. The respective velocity components in these co-ordinates are U and V. The length and velocity 



Effect of transverse curvature on secondary instabilities 



3 




Fig. 2 (a) Velocity profiles compared to those of Tutty et. al. llll at i? = 10"*. The ordinate gives z = ^jR/x*{r — 1). The lowest 
curve is at X = 10^ and the topmost curve is at x = 398. The intermediate curves are spaced at intervals of x'-'^. (b) Dimensionless 
boundary layer thickness 8\fR at = 10*. In both figures, the symbols are from 111], while the lines are present results. 




Fig. 3 (a) Velocity profile at different streamwise locations, at a Reynolds number of 4000. The coordinate r* = (r — l)/5. (b) 
Second derivative of streamwise velocity, if'U jdP' . 



scales which have been used for non-dimensionalisation are the body radius, ro, and the freestream velocity, 
t/oo, respectively, so the Reynolds number is 

(4, 

The solution is obtained by a 3-level implicit finite difference scheme on a uniform grid. At the leading edge, 
two levels of initial data are provided, and downstream marching is commenced. The discretised equation is 
second order accurate in Ax and Ar, and is unconditionally stable in the von Neumann sense. A fairly fine 
grid in the r direction is necessary to capture the velocity and its derivatives accurately. With a grid size of 
10~^ in the x and r directions the results are accurate up to 7 decimal places. 

Velocity profiles at a Reynolds number of 10000 are seen in figure [2| a) to be in good agreement with the 
results of Tutty et. al.[ll|. The dimensionless boundary layer thickness 5 (= ro.99 — 1, where Ur^^^ = 0.99) 
is plotted in figure |2fb). When scaled by the local boundary layer thickness, there is not much of a difference 
visible in the profiles, as seen in figure [3f a) where the Reynolds number is 4000. Here the coordinate r* 
is measured from the body surface. A marked difference near the wall is however evident in the second 
derivative of the velocity (figure|3jb)). This difference is seen below to significantly affect stability behaviour. 
Specifically, an increasingly negative second derivative is indicative of a fuller, and therefore more stable, 
profile downstream. 

The boundary layer does not obey similarity, since there are two parameters in the problem, R and the 
surface curvature defined below. Defining 



S,= {AvxalUy^)^l^, 



(5) 
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where the subscript d denotes a dimensional quantity, we may convert the partial differential equation (HJ to 
an ordinary differential equation in the variable X = r^ (r is non-dimensionalised with body radius ro.) 

l/' + (l + ^^)/=0. (6) 

Here ^ = 2U, and the explicit dependence on x is contained in Sx- It is evident that the velocity profiles would 
be self-similar were the quantity x/R to be held constant. The momentum thickness in an axisymmetric 
boundary layer is of the form 



e = ~ra + Jrl + 2I, where 1= ( U{\-U)rddrd. 

* Jrn 



(7) 



The displacement thickness may be similarly defined. The surface curvature, i.e., the ratio of a typical bound- 
ary layer thickness to the body radius is defined here as 

So = -- (8) 



3 Linear stability analysis 

Based on present wisdom, and our own experience in boundary layer flows, we make the assumption that non- 
parallel effects are small. The equations in this section are the same as those of Tutty et. al. 11111 . expressed 
in terms of the variables introduced by Rao (1967). Flow quantities are decomposed into their mean and a 
fluctuating part, e.g. 

\tot = U{r)x + y{x,r,y,t) (9) 

where v = mx + vr + wy, y being the azimuthal coordinate. Disturbance velocities are expressed in terms of 
generalized stream-functions i// and (/) as 

1 1 fd\if d4>\ ^ d<j> 

u=-^,v = 3^ + 3^ and w = 3^^. (10) 

r or r \ dx ay J or 

In normal mode form 

((//,(/>) (jc,r,y) = ^(^(f,*)(r)exp[i(ax + ny-ft)0]+c.c). (II) 

Here ^{r) and f (r) are the amplitudes of the disturbance stream-functions, a is the wave number in the 
streamwise direction and n is the number of waves encircling the cylinder. The value of n is positive or 
negative for anti-clockwise or clockwise wave propagation respectively. In the temporal stability analysis 
carried out here, the imaginary part of the frequency (0 gives the growth rate of the disturbance. 

Linearising the Navier-Stokes for small disturbances and eliminating the disturbance pressure results in 
two fourth-order ordinary differential equations in f and <^ , given by 



{U - c) (f " - SW - a^W) - W{U" - SU') -{U-c) ancp - - {U"4> + U'^' - SU'^) ) 



I 

iaR 



f "' - 25f "' + 35^^" - 35^ f' - 2a^{W" - SW) + a^W 



-n^S^ (f " - 3SW' -a^W)- na[(P" +S(P'-{a^+ n^S^)(P] 
and 



(12) 



{U - c) + S(P' - n^S^-(P) + U'4>' -{U-c) naS^W -- 



I 

iaR 
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-naS^W" + 3anS^W - {4anS'^ - O^nS^ - an^S^)W 

Here S = d/r, and the boundary conditions are 
f = f'=0, 4>^<P'^0 at r=l, 
and 

W = W'^0, 4> = 'P'^0 as r^oo. 



(13) 



(14) 



(15) 



Upon putting 5 = and letting n ^ oo such that nS tends to a finite quantity corresponding to the spanwise 
wavenumber, j3, equations [12] and[l3]reduce with some algebra to the three-dimensional Orr-Sommerfeld and 
Squire's equations for boundary layers on two-dimensional surfaces (see e.g. |9|). 

The rates of production W+{r) and dissipation W-{r) of disturbance kinetic energy are given by 



W+ir) 
and. 



1 



dU 
dr ' 



(16) 



W-ir) = — (a^luu* +VV* +ww*)+u'u'* +v'v'* +w'w'* 
Re V 

+ ^[n^uu* + {\ + n^){vv* + WW*) + 2in{v*w ~ vw*)]^ 

where the superscript * denotes the complex conjugate. Note that the last term in[l7]is derived from squares 
of magnitudes, and is thus real and positive. 



(17) 



3.1 Inviscid stability characteristics 

It is instructive to first study what happens under inviscid conditions. For two-dimensional flow, the existence 
of a point of inflexion in the velocity profile is a necessary condition |8 | for inviscid instability. The axisym- 
metric analog of this criterion has been derived for various situations e.g.. Duck [2| obtained the generalised 
criterion for axisymmetric disturbances on axisymmetric compressible boundary layers. 

In brief, in the inviscid limit we may eliminate all variables except v in the momentum and continuity 
equations for the linear perturbations, to get 



{U-c) 



(3n^ + ah^) v' 



U"- 



U' 



v = 0. (18) 



From a procedure similar to that for two-dimensional flows, a necessary condition for instability, that the 
quantity / = [/" — [{a^r^ — n^) /r / {a^r^ + n'^)]U' , has to change sign somewhere in the domain, is obtained. 
Letting r ^ oo, we recover the two-dimensional Rayleigh criterion. 

Unlike in two-dimensional flows, the quantity / depends on the streamwise and azimuthal wavenumbers, 
but in order to check for instability it is sufficient to evaluate the limiting cases /i and I2 respectively for 
a/n^O and «/a ^ 0. Using equations [T] and |6] h and Ij can be written as 



/, = U" 



U' 



U' 



and l2 = U" + — = - 



ARxd/ro 



88 ■ 



(19) 
(20) 
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At the wall and at the freestream, g and g" are equal to zero, so h is zero too. In between, I2 is always 
negative since both g and are positive, /j is negative everywhere as well, i.e. / never changes sign. In figure 
|4]these quantities are plotted for a sample case (^o = 0.8). We conclude that the incompressible axisymmetric 
laminar boundary layer on a circular cylinder is inviscidly stable to axisymmetric and non-axisymmetric 
disturbances at any curvature. 

In two-dimensional boundary-layers, the inflexion point criterion has provided a general guideline for 
viscous flows as well, since a flow with a fuller velocity profile typically remains stable up to a much higher 
Reynolds number. We may therefore expect from figure |4] that an axisymmetric boundary layer will be more 
stable than a two-dimensional one. Also as the curvature increases (not shown) the tendency to stabilise will 
be higher. Note that a change of sign in / may occur on converging bodies. We do not consider that case here, 
but mention that the axisymmetric analog of Fjortoft's theorem. 



{U-U, 



u" 



ria^r 



2„2 _ 



-u' 



<o, 



(21) 



where Us is the velocity at the inflection point, being a stricter criterion than the Rayleigh could then be used. 
The above may easily be obtained again by a procedure similar to that in two dimensions. 



3.2 Numerical method and validation 

Equations 1 1 2lto[T5]form an eigenvalue problem, which is solved by a Chebyshev spectral collocation method. 
The transformation 



{\+yc{i))K 
1 + ^-t:^ ^— 



where 

yc{i) =cos{-), i = 0,l...,N. (23) 

are the collocation points, is used to obtain a computational domain extending from r=ltor = L+ l and to 
cluster a larger number of grid points close to the wall, by a suitable choice of K. We ensure that L is at least 
5 times the boundary layer thickness, so that the far-field boundary conditions are applicable. Eigenvalues 
obtained using 81 and 161 grid points are identical up to the sixth decimal place. 
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Tutty et. al. [Ill 






Present 




n 




Rc 


a. 


Cr 




Rc 




Cr 





47.0 


12439 


2.730 


0.317 


47.0 


12463 


2.730 


0.318 


1 


543.0 


1060 


0.125 


0.552 


581.0 


1013 


0.115 


0.552 


2 


91.1 


6070 


0.775 


0.442 


91.0 


6093 


0.775 


0.421 


3 


43.4 


10102 


1.600 


0.403 


43.0 


10110 


1.580 


0.410 


4 


26.8 


13735 


2.540 


0.398 


27.0 


13742 


2.520 


0.401 



Table 1 Critical Reynolds number and other parameters for different modes, in comparison with 1 11]. The streamwise location 
where instability first occurs is denoted as x,.. Uc and c,- are the streamwise wavenumber and phase speed corresponding to the 
critical Reynolds nur 




We compare our critical values with those of Tutty et. al. 1111 in table[T] and find them to be in reasonable 
agreement. The helical mode (n = 1) is destabilised first at a Reynolds number of 1013, and x = 581. The 
axisymmetric (n = 0) mode is unstable only above a Reynolds number of 12463. 

It is only the helical (n = 1) mode which is unstable over a significant axial extent of the cylinder. Even 
this mode is never unstable for curvatures above So= 1, as may be seen in figure|6l At curvature levels below 
this, as well as at low Reynolds numbers, the helical mode is expected to decide dynamics, since the ranges 
of instability of other modes are subsets of the range of the n = I mode. 

Figure |7ja) shows the downstream variation of the critical layer location Vcr, where U (yc,-) = Cr. It is 
seen that as one moves downstream, i.e., as the curvature increases for a given Reynolds number, the critical 
layer moves closer to the wall. Since the production layer scales as 2n{ycr + l)''o, as the cylinder becomes 
thinner and thinner, the cross-sectional area over which production is possible is much smaller, explaining 
the stabilization at large curvature. The energy budget shows that the production layers of several unstable 
modes overlap, and this could give rise to interactions between them, so the nonlinear stability could be very 
different from that in a planar boundary layer. 



4 Secondary instabilities 

A laminar flow containing linear disturbances of a significant amplitude is unstable to secondary modes. The 
A -structures seen in |6| and Kachanov |5 i, considered to be the precursors of turbulent spots, are a signature 
of these modes. As a rule of thumb, nonlinearity in boundary layers becomes detectable when the amplitude 
of the linear (primary) disturbance is 1% of the mean flow. 
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Fig. 6 Critical Reynolds number as a function of curvature for the mode n = 1. 
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Fig. 7 Height of the critical layer as a function of the streamwise distance, n = 1, a = 0.125. Dashed line: Re = 2000; solid line 
Re = 10,000. 



The approach we follow is standard, as in fS^. The periodic basic flow is expressed as, 
Vbasic = U(r)+ApVp 



(24) 



where we have introduced a subscript p for the primary (linear) disturbance. Ap is the ratio of the amplitude 
of disturbance velocity to the freestream velocity. Vp is the disturbance velocity of primary modes, obtained 
from linear stability analysis. The secondary disturbance, in normal mode form, is 



Vs = ^ {\+{r)cx\)[i{k+x + mj^Q - (0+t)] + v_(r)exp[i(^_A; + m_0 - (0-t)] +c.cj , 



(25) 



where k+ and +k- are the streamwise wavenumbers of the secondary waves. The azimuthal wavenumbers of 
secondary waves are m+ and . 

Equations for the secondary instability are obtained by substituting [24] and |25] into the Navier-Stokes 
equations, retaining linear terms in the secondary disturbances, and deducting the primary stability equations. 
The streamwise component of the velocity u+ and m_ are eliminated using the continuity equation. The final 
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equations contain v, w and p. On averaging over x, 9 and f, only the resonant modes survive, which are related 
as follows: 



k^ + k_ = a, in^+m^=n and {(0+ + (0^) ^ = CO . 
The final secondary instability equations are 



(26) 



{{U-c){S- 



D)-U'- 



k+R I 



S[kl + (ml- 



I ([/ - c)im+S + {ki + (ml -l)S'^ + SD~D^^^w+- ik+p+ 



-— — k+Sup + iS Vp — mS Wp + k-u'p + [k+Up — iSvp — m-Swp)D — iVpD 



mS 



— In -m_ I Un H v„ H w„ VnZ) 

2 V k- J 2k- P 2k- 2k- P . 



2im+S 



V+ + 



ik+{U ~c) + -{kX + {mX + l)S''-SD- 



HI-') 



. 1 1 / « 
1 '^/'+2^p + ^^/"^. 



1 



ik^ 



iS 



rri-a 



ik+{U -c) + -[k% + {m\ + \)S-SD- D' 



2im+S^ 

1 

7? 



w*_=0 
ik- 



-P+ 



(27) 



(28) 



imS 



as 
2kl 



Vp 1 a \ 1 
Vo + 4 + -(l + ^jvpZ) 



2 2 



iS t am- 



jVp- ^Jw_ =0 



i2r--. 



(29) 



with three corresponding equations in vl, wl and p*L. The operator D stands for differentiation with respect 
to the radial coordinate. 

The boundary conditions are 



Vs = at r = 1 , Vs ^ as r ^ oo, and p ^0 as r 



(30) 



For the flow under consideration, the growth/decay rates of primary modes are small, hence d {A p) /dt can 
be neglected during one period of time. As mentioned above. Squires theorem does not apply and therefore the 
primary modes must be taken to be three-dimensional. Equations|27]to|29]reduce to the secondary instability 
equations of a flat plate boundary layer by letting 5 = 0, m+S j3, in^S —j5 and ns = 0. The system is 
solved as before. Disturbance growth rates for a zero pressure gradient boundary layer agree well with those 
of El. 



4.1 Results 



The main finding is that for high levels of curvature the flow is stable to secondary modes (as well as the 
linear modes), but secondary modes can extend the curvature range over which disturbance growth is possible. 
As in the case of a two-dimensional boundary layer subharmonic modes are dominant here too. The axial 
wavenumber of the least stable secondary mode is k+ = k- = a/2. At other values fc+ or the growth is 
smaller. This is similar to the behaviour in flat plate boundary layers. The most unstable secondary modes are 
of opposite obliqueness, with azimuthal wavenumber = 2n and m_ = — «. 

is presented in figure[8h. The amplitude Ap of the primary wave is taken to be 2% of [/«,, but the answers do 
not depend qualitatively on this choice. The fiow is seen to be unstable to secondary modes under conditions 
where all primary disturbances decay. For comparison the growth rate of the least stable primary disturbance 
(a = 0.125 and « = 1) is shown as a dotted line. At small primary wave numbers the growth rate is found 
to be small but for a wider band of streamwise locations. As the wavenumber increases the growth rate also 
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0.833 1.196 ° 1.473 1.703 1.905 




2000 4000 

X 



Fig. 8 Growth rate of secondary disturbance waves along the cylinder axis for the non-axisymmetric mode n = 1, m_ = — 1 
with Ap = 0.02 axR= 1000. The most unstable primary mode is shown by the dashed line. 




Fig. 9 Growth rate of secondary disturbance along the cylinder axis for (a) R = 2000 and (b) R = 5000. The other parameters 
are the same as in figure [s] 



increases and instability is restricted to a narrow streamwise region. The maximum growth occurs when 
a = 0.30, and k+ = = 0.15. As discussed earlier m+ = 2 and m- = — 1 in this case. It is found that 
the secondary modes also decay at higher curvature, but the extent of curvatures for which the instability is 
sustained is much larger. 

In figure |9] the growth of secondary waves at higher Reynolds numbers are plotted. The flow conditions 
are same as in figure[8]except the Reynolds numbers. In|9^ Reynolds number is 2000 whilejOj) is for Reynolds 
number 5000. The primary modes amplify substantially in these cases. The behaviour at a higher Reynolds 
number, as seen in figure |9] is as expected. 

Incidentally a small growth is found at small curvatures (figure |9j)) at this Reynolds number . However 
the this growth is very smaller compared to that at higher locations. 

The least stable secondary modes for other values of the azimuthal wavenumber n are shown in figure [TOl 
It is clear that in the range of Reynolds numbers of interest, these modes are not expected to dominate. The 
axisymmetric mode is not shown, but does not afford any surprise either. We have computed growth rates for 
different modes ranging n = I to « = 4, but have shown only the results for n = I and n = 2. There is no 
qualitative difference in other two modes, where the relationship of m+ to n holds good. 



5 Conclusions 

The boundary layer in the flow past a cylinder is stable to (linear and) secondary disturbances at curvatures 
higher than ~ 0{l), i.e. when the radius of the body is of the order of or less than the local boundary 
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Fig. 10 Growth rate of secondary disturbance modes with Ap = 0.02. (a) n = 2, ra+ = 4 and m_ = —2 al R = 3000. (b) n = 3, 
m+ = 6 and ra_ = -3 at i? = 5000. 



layer thickness (shown here in terms of momentum thickness). This indicates that a turbulent axisymmetric 
boundary layer, especially over a thin body, could have a tendency to relaminarise downstream. The flow is 
inviscidly stable at any curvature. Transverse curvature thus has an overall stabilising effect, acting via the 
mean flow and directly through the stability equations. 

The production layers of the disturbance kinetic energy of these modes have a significant overlap, which 
gives rise to the possibility of earlier development of nonlinearities. Thus, while transverse curvature delays 
the first instability, it can contribute once instability sets in to a quicker and different route to turbulence. 

Secondary disturbances remain unstable at larger curvatures than linear modes. However there is again a 
maximum curvature (^o w 2 for Ap = 0.02) above which all disturbances decay. The most unstable secondary 
modes are always those whose azimuthal wavenumbers are related to that of the primary mode by m+ = 2n 
and ifi- = —n. For a helical primary mode (which is the most unstable) this means that one of the secondary 
perturbations is helical as well, but of opposite sense, while the other has the same sense but two waves 
straddle the body. We contrast this to a Blasius boundary layer, where the most unstable secondary mode is 
three-dimensional while the most unstable primary is two-dimensional. There, the spanwise wavenumber j5 
is of the same order as the streamwise wavenumber, and two sets of identical looking waves travel in the 
positive and negative spanwise directions. We do not yet have an explanation for our observation, except to 
say that in a coordinate moving with the primary wave, the observed most unstable secondary is the simplest 
combination containing one forward propagating (in the azimuthal direction) and one backward propagating 
wave relative to the primary. As for the axial wavenumber, the subharmonic modes are least stable, as in 
two-dimensional boundary layers. 

Our studies indicate that experimental and numerical studies of this problem could uncover new physics 
about the transition to turbulence. We hazard a prediction that this process will be significantly different from 
two-dimensional flow. 
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